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The relationship between bats and coronaviruses (CoVs) has received considerable attention since the severe acute respiratory 
syndrome (SARS)-like CoV was identified in the Chinese horseshoe bat (Rhinolophidae) in 2005. Since then, several bats 
throughout the world have been shown to shed CoV sequences, and presumably CoVs, in the feces; however, no bat CoVs have 
been isolated from nature. Moreover, there are very few bat cell lines or reagents available for investigating CoV replication in 
bat cells or for isolating bat CoVs adapted to specific bat species. Here, we show by molecular clock analysis that alphacoronavi- 
rus (a-CoV) sequences derived from the North American tricolored bat (Perimyotis subflavus) are predicted to share common 
ancestry with human CoV (HCoV)-NL63, with the most recent common ancestor between these viruses occurring approxi¬ 
mately 563 to 822 years ago. Further, we developed immortalized bat cell lines from the lungs of this bat species to determine if 
these cells were capable of supporting infection with HCoVs. While SARS-CoV, mouse-adapted SARS-CoV (MA15), and chime¬ 
ric SARS-CoVs bearing the spike genes of early human strains replicated inefficiently, HCoV-NL63 replicated for multiple pas¬ 
sages in the immortalized lung cells from this bat species. These observations support the hypothesis that human CoVs are capa¬ 
ble of establishing zoonotic-reverse zoonotic transmission cycles that may allow some CoVs to readily circulate and exchange 
genetic material between strains found in bats and other mammals, including humans. 


B ats are known to be reservoir hosts for several human viruses, 
including rabies, Marburg, Nipah, Hendra, and the severe 
acute respiratory syndrome coronavirus (SARS-CoV) (5). In ad¬ 
dition, virome studies have shown unprecedented numbers of vi¬ 
ruses present in the fecal samples of this ancient mammalian spe¬ 
cies (11, 24). However, little is known about the genetic 
architecture of most bat species, the virus variation and gene flow 
that occur through different species, the potential of different bat 
species to support human virus replication, the differences be¬ 
tween the bat and human immune systems, or the potential of bat 
viruses to undergo zoonotic transmission to humans and other 
mammals. 

Coronaviruses are the largest known RNA viruses; these 
viruses contain single-stranded plus sense genomes and are clas¬ 
sified in the family Coronaviridae, which is divided into three gen¬ 
era, including Alphacoronavirus (a-CoV), Betacoronavirus (p- 
CoV), and Gammacoronavirus (y-CoV). Five CoVs are known to 
cause human disease, including the P-CoVs SARS-CoV, human 
CoV (HCoV)-OC43, and HCoV-HKUl and the a-CoVs HCoV- 
229E and HCoV-NL63 (35). Three of these HCoVs have been 
shown to or have been predicted to have spilled over from zoo¬ 
notic reservoirs, including SARS-CoV, which likely emerged from 
the Chinese horseshoe bat (Rhinolophidae) (26), HCoV-OC43, 
which probably emerged from bovine CoV (BCoV) (50), and 
HCoV-229E (36), which was predicted by molecular clock analy¬ 
sis to share a most recent common ancestor (MRCA) just over 200 
hundred years ago with a bat CoV found in the leaf-nosed bat 
(Hipposideros caffer ruber) bat in Ghana (36). The close link be¬ 
tween bat and human CoVs has led to the speculation that all 
human, and perhaps mammalian, CoVs may have originated in 
bats (19, 36, 49). 


HCoV-NL63 was first discovered in 2004 as a new HCoV iso¬ 
lated from a 7-month-old baby suffering from bronchiolitis (48). 
A similar virus was isolated around the same time in samples de¬ 
rived from an 8-month-old baby with pneumonia (16). Since 
2004, this HCoV has been detected in 1.0 to 9.3% of respiratory 
tract samples collected from several different countries (15), 
which indicates that HCoV-NL63 is distributed worldwide. There 
is no known reservoir for this virus, and little is known about its 
evolutionary history prior to 2004, although phylogenetic evi¬ 
dence suggests that HCoV-NL63 has infected humans for centu¬ 
ries, as it is predicted to have diverged from HCoV-229E approx¬ 
imately 1,000 years ago (38). 

There are more than 1,100 species ofbats, with bat populations 
inhabiting every continent except Antarctica. Bats belong to the 
order Chiroptera and are further divided into the suborders of 
Yinpterochiroptera, which include mostly the megabats (formerly 
known as the megachiroptera) and Yangochiroptera, which in¬ 
cludes most of the microbat (formerly known as the microchirop- 
tera) families (46). In the United States, most bat species are small 
nocturnal bats that belong to the family Vespertilionidae, and all 
are insectivores. In Maryland, previous surveys have yielded fecal 
samples from the big brown bat (Eptesicus fuscus) , the little brown 
myotis ( Myotis lucifugus), the northern long-eared myotis ( Myotis 
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septentrionalis ), the eastern small-footed myotis ( Myotis leibii), 
the hoary bat ( Lasiurus cinereus), the red bat ( Lasiurus borealis ), 
and the tricolored bat ( Perimyotis subflavus). In our previous 
work, we have shown that five of these species shed unique a-CoV 
sequences (11); however, the reagents necessary to isolate these 
viruses and determine if North American bat CoVs pose a threat 
to public health are currently lacking. 

In this study, we identified nucleic acid sequences that poten¬ 
tially indicate the presence of a novel a-CoV that is predicted to 
share an RCA with HCoV-NL63 in the tricolored bat. We devel¬ 
oped an immortalized lung cell line for this U.S. bat species and 
show that human CoVs grow in these bat cells. This observation 
suggests that human CoVs are capable of infecting multiple mam¬ 
malian hosts, potentially establishing zoonotic-reverse zoonotic 
cycles that allow CoVs to maintain viral populations in multiple 
hosts, evolve novel recombinant viruses with viral genes derived 
from human and animal CoVs, and traffic back into human pop¬ 
ulations at a later date. In addition, these results suggest that 
HCoV-NL63 may have originated in bats and crossed the species 
barrier to infect humans roughly 563 to 822 years ago. Supporting 
a growing body of literature, our data support the hypothesis that 
cross-species transmission events and the emergence and coloni¬ 
zation of new species represent common features of the Corona- 
viridae. 

MATERIALS AND METHODS 

Bat sequences for determining bat phylogeny. All of the available cyto¬ 
chrome b genes ( 1 ) from the North American bats found in Maryland and 
other bats of interest were downloaded from GenBank along with the 
same genes from several other common mammals. The nucleotide gene 
sequences were then aligned by ClustalX, a maximum likelihood tree was 
generated using PhyML with 100 bootstraps, and the tree image was ed¬ 
ited and exported using the bioinformatics tools available in the Geneious 
software suite version 5.4.3 (13). 

Identification of novel a-CoVs in North American bats. In our pre¬ 
vious studies, we demonstrated that a-CoV sequences are present in the 
fecal samples of eastern North American bat species (11). Using the exact 
procedure as previously described ( 11), we used Roche 454 sequencing to 
determine the viral sequences present in bat fecal samples from big brown 
bats captured in the Saratoga National Historical park in New York (New 
England CoV [NECoV]) and tricolored bats from the Chesapeake and 
Ohio Canal National Historical Park in Maryland (Appalachian Ridge 
CoV strain 2 [ARCoV.2]). Then, we used previously reported primers and 
protocols (1 1 ) to amplify a >2,200-nucleotide (nt) fragment in the repli- 
case region of these viruses, encompassing a portion of nsp 13, all of nsp 14, 
and a portion of nsp 15. The amplified fragments were electrophoresed on 
a 1% agarose gel, and the >2,200-nt band was excised, purified, and 
subjected to Sanger sequencing as previously described (11). These se¬ 
quences were deposited into GenBank (see below). 

Phylogenetic and molecular clock analyses of a-CoVs found in 
North American bats, (i) Phylogenetic analysis. The sequences of the 
>2,200-nt fragments of ARCoV.l, ARCoV.2, and NECoV were com¬ 
pared to the same region of several known CoV sequences downloaded 
from GenBank. The sequences were aligned using ClustalX as imple¬ 
mented in Geneious 5.4.3 (13), and the alignment was manually trimmed 
and corrected to generate a 2,321-nucleotide alignment. A maximum like¬ 
lihood tree was generated using PhyML with 100 bootstraps, and the tree 
image was edited and exported using the bioinformatics tools available in 
the Geneious software suite version 5.4.3 (13). This was the largest frag¬ 
ment available for all three genomes, and that was the basis for generating 
the tree using these sequences. 

(ii) Molecular clock analysis. Molecular clock analysis was conducted 
using BEAST version 1.7.1 (14), following the same protocol as that used 


by Pfefferle et al. (2009) (36) and using the same 650- to 800-nt fragment 
of the replicase region of several known CoVs to estimate the date of the 
most recent common ancestor for ARCoV.l and ARCoV.2. The replicase 
sequences for ARCoV.l (11), NECoV, and ARCoV.2 were derived from 
sequence reads obtained by 454 sequencing, and because NECoV and 
ARCoV. 1 were nearly identical, only the ARCoV. 1 sequence was used in 
the analysis. Of note, this sequence is a portion of the viral replicase gene 
(nspl2), which is arguably the most conserved region of the CoV genome, 
making it the most appropriate target for molecular clock analysis. These 
replicase fragment sequences were deposited in GenBank (see below). 
Most of the sequences were dated in years before present, which was 2011 
when this study was conducted. Using the date found by Vijgen et al. 
(2005) (50) for the HCoV-OC43 and bovine CoV sequences, and follow¬ 
ing the method of Pfefferle et al. (2009) (36), a normal probabilistic prior 
with a mean of 121 years before the present time and a standard deviation 
of 13 years was used to calibrate the analysis (36, 50). Both the 
GTR+Gamma 4 + 1 and the SRD06 models were tested under the as¬ 
sumption of an uncorrelated lognormal clock and a constant population 
size. In addition, an exponential population size assumption was tested 
using the SRD06 model. Markov chain Monte Carlo chains were set to 
200,000,000 iterations with sampling every 20,000 generations. Bayes fac¬ 
tor analysis, which was used to compare the models, was conducted using 
the program Tracer version 1.5 (http://tree.bio.ed.ac.uk/software/tracer). 
The SRD06 model was found to be superior to the GTR model by a log 10 
Bayes factor of 240. 

Bat tissue samples. Bats for this study were obtained from the Save 
Lucy Campaign, an education and bat rehabilitation center in Virginia. 
This organization specializes in rehabilitating bats that are sick or injured, 
and all of the bats in the rescue were observed for symptoms of rabies virus 
on intake. Only nonrabid bats that were badly injured beyond recovery 
were euthanized for this study. These included two big brown bats (EpFu) 
and one tricolored bat (PeSu). The bats were euthanized at the rescue 
center, and the bat organs were surgically removed and placed into pro¬ 
cessing medium (calcium- and magnesium-free PBS supplemented with 
100 mg disodium EDTA, 1% penicillin/streptomycin, and 1% gentami¬ 
cin) (7). Organs harvested from each bat included heart, lung, intestine, 
brain, kidney, liver, and stomach. The organs in medium were delivered 
on ice within 8 h to our laboratory for processing. 

Biosafety considerations. All of the bat organs harvested at the bat 
rescue were transported to our lab at the University of North Carolina 
(UNC). The tissues were further processed in a separate and locked bio¬ 
safety level2+ (BSL2+) room that required all personnel to have received 
a prophylactic series of rabies vaccines. All work was performed under 
certified biosafety cabinets, and all personnel wore gloves and Tyvek 
aprons (Fisher Scientific) when working with the fresh tissue. 

Generation of primary bat cell lines from North American bats. The 
organs were removed from transport buffer, washed with cold processing 
media, placed in specimen tubes with fresh processing media, and then 
placed on ice. To prepare the tissues, we followed a published protocol by 
Crameri et al. (2009) ( 7) . The tissues from each organ were finely dissected 
using a sterile scalpel and washed with cold processing media. The dis¬ 
sected tissues were then transferred to a 50-ml conical tube, covered with 
cold 0.25% trypsin (Gibco), and incubated at 4°C overnight. The follow¬ 
ing day, the samples were centrifuged at 37°C on a benchtop shaker set at 
200 rpm for 1 h. Supernatants were filtered through a 40-p.m cell strainer 
(Fisher) into a 50-ml conical tube containing 10 ml of fetal calf serum 
(FCS) (HiClone). The cell strainer was then rinsed with phosphate-buff¬ 
ered saline (PBS) to ensure that all of the individual cells were delivered 
through the strainer. Next, the larger pieces of dissected tissues were in¬ 
cubated with 0.25% trypsin (Gibco) at 37°C for 30 min on the shaker, and 
then the supernatants were add to the FCS in the 50-ml conical tube. This 
was repeated until most of the cells were removed from the tissue scaf¬ 
folds. Next, the cells were pelleted at 800 X g for 5 min in a Beckman 
benchtop centrifuge. The cells were then resuspended in primary cell cul¬ 
ture medium (Dulbecco’s modified Eagle medium [DMEM]/F12-Ham’s 
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media supplemented with 15% Fetal Clone II, 1% penicillin/streptomy¬ 
cin, 1% nonessential amino acids, and 1% gentamicin), transferred to T25 
flasks, and incubated at 37°C with 5% C0 2 . Robust cultures of primary 
heart, lung, brain, and kidney cells were established. The cells were fed by 
removing the medium and replacing it with fresh primary medium with 
15% FCS every 3 to 5 days until the cells reached confluence. The primary 
cells were then maintained for 2 to 5 passages to generate cell stocks that 
were cryopreserved and stored at — 140°C. 

Generation of immortalized bat cell lines from North American 
bats. Primary cells that propagated well enough to reach sufficient 
numbers in up to three passages were selected as candidates for im¬ 
mortalization. Cells were prepared in 6-well plates and grown to 70 to 
80% confluence at 37°C with 5% C0 2 in primary cell culture media. 
The cells were then infected, as previously described (18), with lenti- 
virus vectors expressing the human telomerase reverse transcriptase 
(hTERT) gene and the murine Bmi-1 gene (18), an oncogene that 
promotes efficient self-renewing cell divisions. A mixture of hTERT 
and Bmi-1 was created by placing 1.5 ml of each vector into a 15-ml 
conical tube at vector concentrations of >10 6 infectious units/ml. 
Next, we added 6 p.1 of Polybrene at a concentration of 8 p-g/ml to 
increase the infectivity of the lentivirus vectors, and the mixture was 
vortexed gently and stored on ice. Next, we removed the medium from 
the cells in the 6-well plate and washed the cells twice with PBS 
(Gibco). For the cells to be immortalized, we added 1 ml of the 
hTERT-Bmi-l-Polybrene mixture and incubated it for 3 h at 37°C 
with 5% C0 2 . For the control cells, 1 ml of medium was added prior to 
incubation under the same conditions. After the incubation, the su¬ 
pernatants were removed from the cells and discarded into a 10% 
bleach solution. Primary cell culture medium was then added to all of 
the wells, and the cultures were incubated overnight at 37°C with 5% 
C0 2 . On the following day, the cells were transfected a second time 
using the same protocol. The cells were then monitored daily to look 
for morphological changes in the cells and to determine if the immor¬ 
talization process was successful. Passages were conducted when the 
cells reached confluence, which varied by cell type, with the lung cells 
requiring 3 to 4 days between passages. Immortalization was deter¬ 
mined to be successful in those cells that continued to thrive beyond 
the number of passages required for the primary cells to reach senes¬ 
cence and eventual cell death. To determine the extent of the immor¬ 
talization, the infected cells were propagated for several additional 
passages, with the immortalized lung cells from the tricolored bat 
showing the best postsenescence growth phenotype. Several stocks of 
primary and immortalized cells were cryopreserved at each passage 
and then stored at — 140°C for future use. 

Infection of PESU-B5L cells with HCoV-NL63 and SARS-CoV. Im¬ 
mortalized bat lung cells from the tricolored bat (Perimyotis subflavus) were 
seeded onto a 6-well plate and grown in primary cell culture medium at 37°C 
with 5% C0 2 until they reached 70 to 80% confluence and a concentration of 
~3 X 10 5 cells per well. The supernatants were then removed, and the cells 
were infected with 1 ml of wild-type HCoV-NL63 (wtNL63) or HCoV-NL63 
expressing green fluorescent protein (GFP) (NL63gfp) ( 12) or PBS (mock) at 
a multiplicity of infection (MOI) of 1. For the adsorption period, the cells 
were incubated at 32°C with 5% C0 2 for 1 h with agitation every 15 min. At 
the end of the adsorption period, an additional 2 ml of culture medium was 
added to each well, and the infected cells were maintained at 32°C with 5% 
C0 2 and monitored for 8 days. For SARS-CoVgfp (44) and several related 
strains, the experimental infections were conducted under biosafety level 3 
containment following the same parameters as described above, using 
mouse-adapted SARS-CoV (MA15), a SARS-CoV strain adapted to be lethal 
in mice, GD03, a recombinant chimeric SARS-CoV clone bearing the spike 
sequence of a 2003 human clinical isolate, and HC/SZ/61/03, a recombinant 
chimeric SARS-CoV clone bearing the spike sequence of an early human 
isolate (8, 39, 41, 42, 44). These infected cells were maintained at 37°C with 
5% C0 2 and monitored for 5 days. 

Verification of SARS-CoV and HCoV-NL63 replication in PESU- 


B5L. To verify SARS-CoV and its related strains and HCoV-NL63 repli¬ 
cation in the PESU-B5L lung cells, we used primers designed to amplify 
the subgenomic leader containing mRNA transcripts of each virus by 
reverse transcriptase PCR (RT-PCR). The amplified leader-containing 
cDNAs are specific to each HCoV and would be detected only if the virus 
was undergoing active transcription and replication. For HCoV-NL63, we 
used the primer sets and the exact RT-PCR conditions to amplify the 
subgenomic mRNA encoding the nucleocapsid (N) gene as previously 
described (12). The resulting bands were also excised, purified, and sub¬ 
jected to Sanger sequencing at the UNC Genomic Analysis Center to ver¬ 
ify the origins and sequence of the leader-containing cDNAs. For SARS- 
CoV and its related strains, we used primers designed to detect 
subgenomic transcription of multiple viral proteins using primers and 
reaction conditions that have been published elsewhere (53). The result¬ 
ing band for the SARS-CoV membrane gene was excised, purified, and 
sequenced in the same way as described for the HCoV-NL63 N gene. 

Detection of ElCoV-NL63 in PESU-B5L by immunofluorescence as¬ 
say. PESU-B5L cells were infected with HCoV-NL63 at an MOI of 0.7 or 
mock infected. Seventy-two hours postinfection, monolayers were 
washed twice with cold phosphate-buffered saline and fixed with a 50% 
methanol-50% acetone solution. Cells were then permeabilized with 0.1% 
Triton X-100-5% bovine serum albumin (BSA) in PBS and blocked with 
5% BSA in PBS. Rabbit anti-NL63 nucleocapsid sera (the kind gift of Lia 
van der Hoek) were used as the primary antibody. Unbound primary 
antibody was removed by washing three times with 1% BSA-0.05% Non- 
idet P-40 in PBS, followed by addition of an Alexa 546-conjugated goat 
anti-rabbit IgG (Invitrogen). After washing off unbound secondary anti¬ 
body, infected cells were photographed using an inverted microscope with 
a mercury lamp light source and a tetramethyl rhodamine isothiocyanate 
(TRITC) filter. 

Detection of EICoV-NL63 in PESU-B5L cells by Western blotting. 

PESU-B5L cells were infected with HCoV-NL63 at a multiplicity of infec¬ 
tion of 0.7 or mock infected. Seventy-two hours postinfection, cells were 
washed with cold PBS followed by lysis in AV lysis buffer (20 mM Tris- 
HC1 [pH 7.6], 150mMNaCl, 0.5% deoxycholine, 1.0% Nonidet P-40, and 
0.1% sodium dodecyl sulfate [SDS]). The lysates were spun down at 
13,000 rpm for 10 min to pellet the nuclei. Equal volumes of 0.9% SDS and 
10 mM EDTA were added to the postnuclear supernatants (to raise the 
concentration of SDS to 0.5% for virus inactivation), and the mixture was 
incubated at room temperature for 15 min to ensure complete inactiva¬ 
tion of HCoV-NL63. Equal amounts of cell lysates were then loaded onto 
a 4 to 12% Bis-Tris Nupage gel (Invitrogen), and postelectrophoresis, the 
proteins were transferred onto polyvinylidene difluoride (PVDF) mem¬ 
branes using an Xcell blot module (Invitrogen) according to the manu¬ 
facturer’s instructions. Membranes were blocked with 5% nonfat dry milk 
in TBS-T buffer (0.1% Tween 20 in Tris-buffered saline) for 1 h and 
incubated overnight with rabbit anti-NL63 nucleocapsid sera. Unbound 
antibody was removed by washing in TBS-T followed by incubation with 
an antirabbit horseradish peroxidase (HRP) secondary antibody. Proteins 
were visualized using a Pierce Super Signal Western blot chemilumines¬ 
cence kit. For the loading control, a purified mouse monoclonal antibody 
to beta actin (catalog no. 612656; BD Transduction Technologies) was 
used, followed by staining with secondary antimouse HRP antibody (GE 
Amersham Biosciences). 

Nucleotide sequence accession numbers. The sequences determined 
in this study were deposited into GenBank under accession numbers 
JX537911 through JX537914. 

RESULTS 

Evolutionary relatedness of North American bats. The mito¬ 
chondrial cytochrome b genes of several bat and other mamma¬ 
lian species were compared to show the evolutionary relatedness 
of North American bats frequently found in Maryland, including 
the big brown bat (EpFu), the little brown myotis (MyLu), the 
eastern small-footed myotis (MyLe), the northern long-eared 
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FIG 1 Phylogeny of the mitochondrial cytochrome b gene from several mam¬ 
malian species. This maximum likelihood tree shows that North American 
bats of the family Vespertilionidae are more closely related to themselves than 
to other bat or common mammalian species. EpFu, big brown bat; PeSu, 
tricolored bat; MyLu, little brown myotis; MyLe, eastern small-footed myotis; 
MySe, northern long-eared myotis. The scale bar represents nucleotide substi¬ 
tutions. Only nodes with bootstrap support above 70% are labeled. 


myotis (MySe), and the tricolored bat (PeSu) (Fig. 1). The North 
American bat species group in the same cluster with other bat 
species, showing a closer relationship between these species. EpFu 


and PeSu were more closely related to each other than to other 
myotine species found in Maryland (Fig. 1). 

Identification and characterization of novel ol-CoV se¬ 
quences in big brown bats and the tricolored bat. Roche 454 
sequencing was used to identify a-CoV sequences in two addi¬ 
tional bat populations, including a big brown bat population from 
Saratoga National Historical Park in New York (NECoV) and a 
tricolored bat population that was sampled at the Chesapeake and 
Ohio National Historical Park in Maryland. RT-PCR was used to 
amplify a >2,200-nt region of these viruses (comprising a portion 
of nspl3, all of nspl4, and a portion of nspl5), and these were 
sequenced and compared to the same fragment that had been 
amplified previously from a big brown bat population sampled in 
Maryland and was known as ARCoV. 1 ( 11 ). Other regions of these 
genomes have been determined by 454 sequencing, but this was 
the largest fragment available in common between all three se¬ 
quences. Interestingly, the two >2,200-nt sequences derived from 
big brown bats were nearly identical (approximately 97% identical 
at the nucleotide level), suggesting that big brown bats harbor 
similar a-CoVs regardless of geographic location (Fig. 2A and B). 
In contrast, the sequence from ARCoV.2 was approximately 74% 
identical at the nucleotide level, suggesting that different U.S. bat 
species harbor unique a-CoV strains (Fig. 2B). Together, these 
fragments of North American bat CoVs suggest a new clade in the 
Alphacoronavirus genus (Fig. 2B). 

Molecular clock analysis. During the annotation stage, many 
of the CoV-like sequences from tricolored bat samples were most 
closely related, albeit distantly, to HCoV-NL63 (Fig. 2B). There¬ 
fore, molecular clock analysis was conducted to estimate an 
MRCA between ARCoV.2 and HCoV-NL63, using a similar data 
set and the same parameters as those employed by Pfefferle et al. 
(2009) (36). This analysis was performed using sequence frag¬ 
ments obtained from the highly conserved viral polymerase 
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FIG 2 A novel a-CoV in the tricolored bat is closely related to HCoV-NL63. Novel a-CoV sequences have been found in the fecal samples of several North 
American bat species. (A) Schematic showing where the two fragments (used for the trees in panels B and C) occur in the NL63 genome. (B) A >2.2-kb fragment 
of the replicase region (starting at position 16480 in the schematic) was sequenced from three samples, including two from the big brown bat (NECoV and 
ARCoV. 1) and one from the tricolored bat (ARCoV.2). A maximum likelihood tree comparing the nucleotide sequences of these bat CoVs to other known CoVs 
indicated that NECoV and ARCoV. 1 are very closely related while ARCoV.2 is significantly different from NECoV and ARCoV. 1. The three novel a-CoV 
sequences form a novel cluster in the a-CoV group. (C) Molecular clock analysis using a 650- to 800-nt portion of the highly conserved replicase region (starting 
at position 13810 in panel A) predicted that the MRCA of bat CoVs from the Hipposideros caffer ruber bats and HCoV-229E was likely to have existed 212 to 350 
years ago (in agreement with Pfefferle et al., 2009 [36]). The MRCA for HCoV-NL63 and ARCoV.2 was predicted to have existed 563 to 822 years ago. NECOV 
and ARCoV. 1 were identical in this region, so only ARCoV. 1 is shown, and it clustered with other bat a-CoVs. 
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TABLE 1 Models used and highest posterior densities 


Model 

Root (range)" 

Mean substitution rate fo (range) 

GTR + G + I, constant 

SRD06, constant 

SRD06, exponential 

3469(1237-5923) 

5940(2424-9921) 

5450(2143-9358) 

1.629 X 10 -4 (1.034 X 10~ 4 -2.3499 X 10~ 4 ) 
1.76 X 10 -4 (1.073 X 10 -4 —2.570 X 10~ 4 ) 
1.847 X 10 -4 (1.078 X 10~ 4 -2.668 X 10~ 4 ) 

a Values are years BCE (before the common era). 
b Mean substitution rate is the number of substitutions per site per year. 


(nspl2), which is arguably the most conserved region of the CoV 
genome (Fig. 2A). Briefly, both the GTR+Gamma 4 + 1 and 
SRD06 models were tested under the assumption of an uncorre¬ 
lated lognormal clock and constant population size. In addition, 
an exponential population size assumption was tested using the 
SRD06 model. Markov chain Monte Carlo chains were run for 
200,000,000 iterations with sampling every 20,000 generations. 
Bayes factor analysis, used to compare the models, was conducted 
using the program Tracer. The SRD06 model was found to be 
superior to the GTR model by a log 10 Bayes factor of 240. The 
exponential population assumption was found to be superior to 
the constant assumption, but not significantly so (log 10 Bayes fac¬ 
tor of 2). The results of the analyses are shown with highest pos¬ 
terior densities (HPD) in Table 1. Interestingly, the analysis reca¬ 
pitulated the results obtained by Pfefferle et al. (2009) (36) and 
provided an estimated MCRA between ARCoV.2 and HCoV- 
NL63 of approximately 1449 CE (range, 1190 to 1449 CE) (Table 
1 and Fig. 2C). This observation suggests that ARCoV.2 and 
HCoV-NL63 originated from the same ancestor, predicting that a 
potential cross-species transmission event occurred about the 
time that Columbus arrived in North America. Because ARCoV. 1 
and NECoV were nearly identical, only ARCoV. 1 was used in this 
analysis, and this sequence clustered with other a-CoVs of bats 
(Fig. 2C). 

Generating primary cells from sacrificed bats. Bat tissues 
used for this study were obtained from The Save Lucy Campaign, 
which operates as a bat rescue in Virginia. Bats that were too badly 
damaged to be rescued were euthanized, and their organs were 
harvested and transported to UNC, where the tissues were surgi¬ 
cally dissected and trypsinized, and the filtered primary cells were 
plated and grown. Once a sufficient number of cultured primary 
cells were available, these were tested for rabies virus by RT-PCR. 
All of the tested bat tissues were negative for rabies virus (data not 
shown). 

A female tricolored bat that arrived in the rescue in mid-July of 
2011, suffering injuries sustained from being hit by a car, pro¬ 
duced the best results. Brain, heart, liver, lung, intestine, stomach, 
and kidney tissues were obtained from this bat, and primary cells 
were successfully generated for the brain, lung, stomach, and 
heart. The other tissue types grew initially, and either they reached 
senescence before reaching confluence or there were too few cells 
to establish long-term growth. The tricolored bat lung cells grew 
most efficiently, and these cells were passaged three times to gen¬ 
erate a large stock of primary lung cells for this bat (Fig. 3A). 
Stocks of these cells were cryopreserved at — 140°C for future use. 

Generating immortalized lungs cells from the tricolored bat. 
Primary tricolored bat lung cells that had been passaged 3 times 
were immortalized by two concurrent infections with two lentivi- 
rus vectors, one containing hTERT and one containing Bmi-1 
(18). Immortalized and primary control cells were then passaged 
at 3- to 4-day intervals to determine if the immortalization process 


extended the viability of those cells infected with the immortaliza¬ 
tion vectors. By passage 10, the primary cells had reached senes¬ 
cence and had stopped growing completely, while the cells that 
were immortalized continued to grow and divide efficiently (Fig. 
3B and C). The immortalized cells were named PESU-B5L cells, 
and these cells have continued to show prolific growth. At the time 
of writing, we have PESU-B5L cells that have undergone over 30 
passages and over 100 doublings. There were no observable mor¬ 
phological changes in these cells over passage, although early pri¬ 
mary cultures had heterogenous mixtures of cell types. However, 
the immortalization process appears to have selected for a single 
fibroblast-resembling cell population (Fig. 3C). Total RNA was 
extracted from the immortalized cells, and RT-PCR and sequenc¬ 
ing were used to determine that the murine Bmi-1 gene was ex¬ 
pressed in the immortalized cells (data not shown). 

Experimental infection of PESU-B5L lung cells with human 
CoVs. To assess whether or not bat cells could support growth of 
HCoVs, we infected the PESU-B5L cells with HCoV-NL63, a hu¬ 
man a-CoV, and SARS-CoV, a human (3-CoV, using engineered 
versions of these viruses that express GFP during active replication 
(12, 17, 44, 53), so that we could monitor the presence of a pro¬ 
ductive infection by fluorescence microscopy. In the case of 
HCoV-NL63, we detected GFP fluorescence at day 2 postinfection 
and it continued until day 8, when the cells were harvested due to 
overconfluence. While a cytopathic effect was observed in some 
cells by day 5 postinfection (Fig. 4), the cultures remained viable, 
suggesting that many cells were resistant to infection. The progeny 



FIG 3 Immortalization of tricolored bat lung cells, PESU-B5L. Primary tri¬ 
colored bat lung cells were grown in enriched growth medium and passaged 
four times to generate a suitable quantity for immortalization. Primary cells 
were immortalized with lentiviral vectors containing hTERT and the Bmi-1 
proto-oncogene. (A) Primary tricolored bat lung cells. (B) Primary PESU-B5L 
lung cells reached senescence by passage 10. (C) Immortalized cells at passage 
10, including 4 passages as primary cells and 6 passages postimmortalization. 
The immortalized PESU-B5L cells were still going strong after more than 30 
passages (over 100 doublings). 
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FIG 4 Infection of PESU-B5L cells with human coronavirus NL63 expressing 
GFP. PESU-B5L cells were infected with mock or NL63gfp at an MOI of 0.25 
and monitored by bright-held and fluorescence microscopy. (A and B) Mock- 
infected cells by bright-held (A) and fluorescence (B) microscopy. (C through 
F) Cells infected with NL63gfp visualized at day 3 postinfection by bright-held 
microscopy (C) and by fluorescence microscopy (D) and at day 5 postinfection 
by bright-held microscopy (E) and by huorescence microscopy (F). 
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FIG 5 Verification of NL63 replication in PESU-B5L cells. Total RNA was 
harvested from cells on day 5 postinfection and subjected to RT-PCR using 
primers designed to detect leader-containing transcripts of the N gene, which 
would be present only if NL63 replicated. (A) Subgenomic transcripts were 
detected in the PESU-B5L cells infected with both wtNL63 and NL63gfp. (B) 
The sgN band for wtNL63 was excised, and the DNA was purified and se¬ 
quenced to verify that the ElCoV-NL63 leader sequence and the 5' end of the N 
gene were present. L, ladder; N, wt-NL63; G, NL63gfp; M, mock. The leader 
sequence has a beige bar beneath it, and the 5' end of the N gene has a green bar 
beneath it. 


virions harvested from the supernatants of these infections were 
viable upon passage to fresh cells and formed plaques on LLC- 
MK2 cells; however, the titer was low (<10 3 PFU/ml), suggesting 
a potential block in viral egress. Attempts to generate HCoV-NL63 
plaques in the PESU-B5L cells were unsuccessful. 

In the case of SARS-CoVgfp, no GFP expression was detected 
in PESU-B5L cells over 5 days of infection, but GFP fluorescence 
was readily detected in Vero E6 cells infected with the same virus 
as a positive control by 24 h postinfection (data not shown), indi¬ 
cating that the stock of SARS-CoVgfp virus was viable. 

Detection and verification of replication of HCoV-NL63 in 
PESU-B5L cells. To verify that HCoV-NL63 replication had oc¬ 
curred in the PESU-B5L lung cells, we used RT-PCR to amplify 
leader-containing transcripts of the subgenomic N gene. The RT- 
PCR resulted in a ladder of amplicons representing subgenomic 
transcripts that had been generated during the course of efficient 
virus replication (Fig. 5A), and these were sequenced to verify the 
presence of leader-containing transcripts during HCoV-NL63 
replication (Fig. 5B). In addition, infection of PESU-B5L cells by 
HCoV-NL63 was detected by immunofluorescence assay (IFA) 
(Fig. 6A) and Western blotting (Fig. 6B), both directed against the 
nucleocapsid protein. For the IFA, cells expressing HCoV-NL63 
nucleocapsid appeared bright green under fluorescent conditions, 
while uninfected cells were not stained. Western blotting of cell 
lysates using the anti-NL63 N antibody revealed that the infected 
cells strongly expressed a protein of approximately 45 kDa, which 
was consistent with the size of the HCoV-NL63 nucleocapsid pro¬ 
tein. 

Detection of replication of SARS-CoV in PESU-B5L cells. De¬ 
spite the fact that SARS-CoVgfp was not detected by fluorescence 
in the infected cultures of PESU-B5L cells, RT-PCR was con¬ 
ducted on RNA harvested from these cells at 24 and 48 h using 


primers designed to detect the subgenomic GFP transcript, and 
this resulted in the detection of a band in the 24-h culture, indi¬ 
cating that subgenomic transcription was ongoing in these cells 
(Fig. 7). No viruses or replication was detected upon passage of 
supernatants from these cultures. 

Detection of replication of MAI 5 and chimeric SARS-CoV 
strains in PESU-B5L cells. In addition to performing infections 
with SARS-CoV expressing GFP, we also conducted experimental 
infections with MAI 5 and two chimeric SARS-CoV viruses that 
encoded the spike proteins from two early human isolates known 
as GD03 and HC/SZ/61/03. All three of these viruses replicated 
subgenomic RNAs in the PESU-B5L lung cells (Fig. 8), as detected 
by RT-PCR using the total RNA of cells harvested at 48 and 72 h 



FIG 6 Evidence of HCoV-NL63 infection in PESU-B5L cells. PESU-B5L cells 
were infected with HCoV-NL63 at an MOI of 0.8. At 72 h postinfection, cells 
were probed for the nucleocapsid protein of HCoV-NL63 by immunofluores¬ 
cence assay and Western blotting. (A) I, mock-infected cells, fluorescence mi¬ 
croscopy; II, mock-infected cells, bright-field microscopy; III, HCoV-NL63- 
infected cells, fluorescence microscopy; IV, HCoV-NL63-infected cells, 
bright-field microscopy. (B) Western blot of HCoV-NL63 and mock-infected 
PESU-B5L cell lysates with (3-actin as a loading control. Nucleocapsid (N) 
protein appears as a distinct band at ~45 kDa. 
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FIG 7 Evidence of SARS-CoV replication in PESU-B5L cells. PESU-B5L cells 
were infected with SARS-CoV at an MOI of ~ 1, and total RNA was harvested 
from the cells at 24-h intervals and subjected to RT-PCR using primers de¬ 
signed to detect leader-containing transcripts of the subgenomic GFP 
(sgGFP), indicative of SARS-CoV replication. Subgenomic transcripts were 
detected at 24 h postinfection but disappeared by 48 h postinfection, suggest¬ 
ing that SARS-CoV replicates in these cells at a low level. L, Promega 1-kb 
ladder; M, mock-infected cells; V, SARS-CoV-infected cells. 


postinfection and primers designed to detect subgenomic tran¬ 
scription of viral proteins. Subgenomic transcripts were detected 
at both time points for all three strains (Fig. 8A). The band repre¬ 
senting the M gene was excised and sequenced by Sanger sequenc¬ 
ing for all three of these strains, and all contained the appropriate 
sequence to verify the presence of the leader sequence combined 
with the M gene (Fig. 8B and data not shown). Interestingly, de¬ 
spite the fact that the PESU-B5L cells infected with MA15 and the 
two chimeric viruses showed moderate amounts of cytopathology 
that occurred prior to cell senescence (data not shown), we were 
unable to detect any signs of viral replication when the superna¬ 
tants from these infections were passed onto fresh PESU-B5L cells. 

DISCUSSION 

Bats are important reservoir hosts for a large number of emerging 
viruses that cause human disease. However, only a small number 
of species (~36 of 1,100) (5, 11,24) have been studied and shown 
to harbor RNA viruses, and many of these viruses share sequence 
similarity to human strains. The fact that CoVs have been found in 
multiple Old World and New World bat species suggests that bats 
may be particularly suited for maintaining CoVs in nature. In¬ 
cluding the findings in this study, at least three human CoVs have 
been linked to bat CoVs (26, 36), and there is mounting evidence 
that all mammalian CoVs may have originated from bats (19, 36, 
49), particularly human coronaviruses. 

Although the exact number of infections is unknown, HCoV- 
NL63 has been shown to be distributed worldwide; it infects a 
significant number of children and elderly persons each year, and 
HCoV-OC43 and -229E account for nearly one-third of all colds. 
Moreover, SARS-CoV infected over 8,000 people worldwide, with 
a mortality rate approaching 10% (34). These observations indi¬ 
cate that CoVs are important human pathogens that cause signif¬ 
icant morbidity and frequent mortality. The fact that over 1,000 
bat species have not been assessed for the presence of CoVs or 
other RNA and DNA viruses may have enormous public health 
consequences in an outbreak setting. Given that there are more 
than 1,100 species of bats and nearly one-half of the 36 bat species 
that have been studied for viruses have shown evidence of CoV 
sequences, there may be hundreds of novel coronaviruses in bats 
throughout the world. If this is the case, then the next CoV spill¬ 
over from bats to humans is just a matter of ecological opportu¬ 
nity. 
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FIG 8 Evidence of SARS-CoV replication in PESU-B5L cells. PESU-B5L cells 
were infected with the mouse-adapted strain of SARS-CoV, MA15, and with 
chimeric SARS-CoV viruses bearing the spike genes of GD03 and HC/SZ/ 
61/03 at an MOI of ~1. Total RNA was harvested from the cells at 48 and 72 h 
postinfection and subjected to RT-PCR using primers designed to detect lead¬ 
er-containing subgenomic transcripts indicative of replication. (A) Sub¬ 
genomic (sg) transcripts for XI, envelope (E), membrane (M), and open read¬ 
ing frame 6 [ORF6 (X3) ] were detected at both time points for all three strains. 
Lane 1, MA15 at 48 h; lane 2, MA15 at 72 h; lane 3, Invitrogen 1-kb Plus ladder; 
lane 4, SARS+GD03 at 48 h; lane 5, SARS+GD03 at 72 h; lane 6, mock- 
infected cells; lane 7, SARS+HC/SZ/61/03 at 48 h; and lane 8, SARS+HC/SZ/ 
61/03 at 72 h. (B) The sg M bands for all three strains were excised, and the 
DNA was purified and sequenced to verify that the SARS-CoV leader sequence 
and the 5' end of the M gene were present. Shown here is sgMfrom SARS-CoV 
strain GD03. 


North American bats have been shown to harbor a-CoV se¬ 
quences (10, 11, 24, 31, 33); however, no a-CoVs have been iso¬ 
lated from these bats, and in fact, no full-length CoV sequences 
have been obtained or reported from North American bats. Of 
note, in our metagenomics study of North American bats in the 
Eastern United States, we identified a-CoV sequences in multiple 
bat species, including the big brown bat and the tricolored bat 
(11). In this study, we detected a-CoV sequences in a big brown 
bat population in New York and from a tricolored bat population 
sampled in Maryland (Fig. 2). Interestingly, the two sequences 
derived from big brown bats (ARCoV.l and NECoV) were nearly 
identical in a >2,200-nt fragment of the replicase gene, even 
though the two populations were from different regions (New 
York and Maryland) (Fig. 2B). In contrast, the same sequence 
fragment obtained from tricolored bats was significantly different 
(74% identity at the nucleotide level), suggesting that different bat 
species harbor CoVs that have adapted specifically to that species 
(Fig. 2B). Whether or not a compartmentalized CoV from one 
species is more likely to infect humans, as occurs with the rabies 
virus (6, 32, 40), is yet to be determined. However, it is clear that 
specific bat species harbor bat SARS-like viruses that are more 
closely related to SARS-CoV (26). 

The fact that some of the Roche 454 sequences derived from 
the tricolored bat were more closely related to HCoV-NL63 than 
any other a-CoV suggested that this novel North American 
a-CoV may be related to HCoVs. Therefore, we conducted mo¬ 
lecular clock analysis using BEAST to determine the relatedness of 
this a-CoV to other CoVs, by adding it to previous data sets used 
to determine ancestry for other CoVs (36, 50). Of note, ARCoV.2 
and HCoV-NL63 appear to have an MRCA that occurred just over 
560 years ago, based on a 650- to 800-nt portion of the highly 
conserved replicase gene. This analysis was limited to this frag¬ 
ment so that we could use the data obtained from previous studies 
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to calibrate the clock and validate our results. If these predictions 
are correct, this observation suggests that HCoV-NL63 may have 
originated from bats between 1190 and 1449 CE. However, it is 
possible that recombination between bat and human CoVs or 
other mammalian strains may have occurred in this region, and 
therefore, more sequence information is necessary to work out the 
evolutionary history of these CoVs. Comparing the two trees in 
Fig. 2 suggests that recombination has likely occurred, as the tree 
topologies are slightly different even in the nonstructural gene 
sequences that are known to be highly conserved among CoVs. 

Immortalized cell lines are important reagents for investigating 
the similarities and differences between viruses that grow in bats 
and other mammals. To date, only a few immortalized bat cell 
lines have been developed, and the majority of these are for fruit 
bats, which allows for further studies with the Nipah and Hendra 
viruses (4, 7, 22). In this study, we successfully immortalized bat 
cells (Fig. 3) from a bat species known to harbor a unique a-CoV 
that appears to share common ancestry with HCoV-NL63. The 
facts that SARS-CoV and related strains replicate (Fig. 7 and 8) 
and that HCoV-NL63 grows in the PESU-B5L immortalized lung 
cell line indicate the importance of developing robust reagents to 
a number of different bat species. To our knowledge, this is the 
first bat cell line that supports growth of HCoVs. We are currently 
in the process of characterizing the innate immune response in 
these cells, and we have primary cells available for three additional 
New World bat species that will be targeted for immortalization. 
These cell lines will be used to generate reagents for directly study¬ 
ing CoVs in bats and for assessing the potential of these viruses to 
infect humans. Our data suggest that the lentivirus-based immor¬ 
talization system reported herein would also be successful for 
multiple North American bat species and likely represents a ro¬ 
bust tool for rapidly establishing continuous cell lines from bats 
and other zoonotic species. 

Only a small number of immortalized bat cell lines currently 
exist, and these were derived both from Old World bats, including 
the Egyptian fruit bat (Rousettus aegyptiacus) (22), the straw-col¬ 
ored fruit bat ( Eidolon helvum) (4), and the black flying fox ( Ptero- 
pus alecto) (7), and from New World bats, such as the Brazilian 
free-tailed bat ( Tadarida brasiliensis) (ATCC CCL-88). In these 
cases, the cell lines were immortalized using three different meth¬ 
ods, including the simian virus 40 (SV40) large tumor antigen (4, 
7), the adenovirus serotype 5 E1A and E1B genes driven from the 
promoters for human phosphoglycerate kinase and thymidine ki¬ 
nase of herpes simplex virus (22), and hTERT (7). Initially, we 
attempted to immortalize the various primary bat cells using 
hTERT, but these attempts did not result in immortalized cells 
(data not shown). Interestingly, it has been shown that hTERT 
alone cannot immortalize primary human bronchial epithelial 
(hBE) cells (30), possibly due to suboptimal culture conditions. 
Moreover, while hTERT in combination with viral oncogenes has 
been used to immortalize mammalian cell lines, these immortal¬ 
izations have often resulted in immortalized cell lines that were 
limited in differentiation capacity and that frequently exhibited 
other morphological and genetic abnormalities (18). To circum¬ 
vent these limitations, we used a viral-oncogene-independent ap¬ 
proach that was successfully used to immortalize hBE cells for the 
study of cystic fibrosis (18). The primary difference in this ap¬ 
proach was that it used the expression of the murine Bmi-1 gene 
(18), which is a proto-oncogene that normally maintains stem cell 
populations but has also been shown to recapitulate normal cell 


structure and function (18). The PESU-B5L cells immortalized by 
this approach have shown no morphological modifications over 
time and appear to be a homogenous population of lung epithelial 
cells. 

The isolation and growth of novel a-CoVs from the North 
American and other global bat populations continue to be a diffi¬ 
cult problem. We have attempted to isolate a-CoVs of several bat 
species from fecal samples that were positive by PCR for an a-CoV 
by inoculating freshly filtered fecal supernatants from these sam¬ 
ples onto a variety of cell types, including PESU-B5L lung cells. To 
date, we have been unsuccessful in isolating a CoV using this ap¬ 
proach. Therefore, we have been using high-throughput and 
Sanger sequencing approaches to identify and fill in genomic se¬ 
quences that can be engineered into infectious clones that will 
allow us to study these viruses in the laboratory. We previously 
used this synthetic approach to resurrect Bat CoV HKU3, the 
SARS-like CoV identified in the Chinese horseshoe bat (26), and 
showed that the only block to replication of this bat virus in pri¬ 
mate cells was the 180-amino-acid portion of the spike glycopro¬ 
tein known as the receptor-binding domain (RBD) (3). 

In this paper, we report the first evidence that traditional hu¬ 
man CoVs are capable of growth in bat cells. Interestingly, both 
HCoV-NL63 and SARS-CoV and its related viruses replicated in 
these cells; however, replication occurred over multiple passages 
with HCoV-NL63. Recent studies have demonstrated that SARS- 
CoV and HCoV-NL63 use angiotensin-1 converting enzyme-2 
(ACE2) as a receptor for entry into permissive cells (25,28,37,45). 
However, crystallography studies have shown that these HCoVs 
interact with ACE2 using different receptor binding domains in 
the spike glycoproteins of each virus to interact with different 
regions of the ACE2 molecule (23, 27, 52). In fact, we would pre¬ 
dict that it is differences in these interactions that impact the abil¬ 
ity of SARS-CoV to efficiently replicate upon passage in these cells. 
Interestingly, while SARS-CoVgfp replicated leader-containing 
transcripts at 24 h, transcription was not detected at 48 h and no 
GFP fluorescence was observed in these cultures. However, MAI 5 
and the two chimeric human SARS-CoVs bearing spike genes 
from the early stage of the SARS epidemic appeared to replicate 
leader-containing transcripts more efficiently and exhibited more 
cytopathology, indicating that these viruses grew better in these 
cells than did SARS-CoVgfp. This observation suggests that the 
spike protein may be the primary determinant for infectivity in 
these cells and that spike proteins that are more adapted to human 
ACE2, such as SARS-CoV, are less capable of establishing infec¬ 
tion in these bat cells. However, HCoV-NL63 appeared to repli¬ 
cate more efficiently than any of the SARS-CoV variants, suggest¬ 
ing that the HCoV-NL63 interaction site on the tricolored bat 
ACE2 may be the primary determinant of infection. The interac¬ 
tions between PESU-B5L cells and various spike proteins are cur¬ 
rently being explored in more detail. Interestingly, a recent study 
has shown that signatures of recurrent positive selection in the bat 
ACE2 gene map almost perfectly to known SARS-CoV interaction 
surfaces, suggesting that ACE2 utilization preceded the emergence 
of SARS-CoV-like viruses from bats (9). 

In addition, the sequence of the orthologous ACE2 receptor of 
the tricolor bat has not been determined but likely presents the 
determinants of cross-species transmission, allowing HCoVs to 
infect bats. Supporting this hypothesis, several bat ACE2 receptors 
were recently shown to function as receptors for SARS-CoV dock¬ 
ing and entry, and notably, HCoV-NL63 receptor usage was not 
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evaluated in these studies (21). Of note, most a-CoVs use amino- 
peptidase N (APN) molecules as the cellular receptor, and the 
exogenous expression of feline APN in cells that are refractory to 
infection has rendered these cells susceptible to most a-CoVs 
(47). The fact the HCoV-NL63, which is an a-CoV, uses ACE2 is 
an unusual feature of this HCoV, although it is not know which 
receptors the bat a-CoVs use. We are currently in the process of 
cloning and sequencing all of the known CoV receptor ortho- 
logues for this bat species to determine which receptor(s) is used 
by SARS-CoV and HCoV-NL63 and to assess the risk of reverse 
zoonotic transmission from humans to bats. We also plan to eval¬ 
uate PESU-B5L infectivity and receptor usage by other coronavi- 
ruses such as HCoV-229E and EIKU3 in future experiments. At¬ 
tempts to obtain these viruses for this study were unsuccessful. 

CoVs have a long history of cross-species transmission (35). 
BCoV and HCoV-OC43 are closely related, with an MRCA pre¬ 
dicted to have occurred ~100 years ago (2, 50, 51). BCoV strains 
have also spread to alpaca and wild ruminants. These observations 
suggest that HCoV-OC43 arose from a cross-species transmission 
of BCoV into humans, although it is equally likely that the trans¬ 
mission happened in reverse (2, 50, 51). In the a-CoV genus, 
canine CoV (CCoV), feline CoV (FCoV), and the porcine CoV 
transmissible gastroenteritis virus (TGEV) carry genetic evidence 
of recombination with each other, suggesting that these CoVs 
likely originated from or infected the same host prior to becoming 
host range restricted. In fact, early strains CCoV-1 and FCoV-1 are 
thought to have diverged from a common ancestor, and multiple 
recombination events between these strains and an unknown CoV 
(in an unknown host) likely resulted in the novel CoVs CCoV-II 
and FCoV-II (20a). Sequence analysis looking at the similarities 
between CCoV-II and TGEV suggests that TGEV emerged from a 
cross-species transmission of CCoV-II from an infected canine to 
pigs (29). 

SARS-CoV emerged from the Chinese wet markets, where the 
virus is thought to have crossed the species barriers from Chinese 
horseshoe bats to masked palm civets (Paguma larvata) and rac¬ 
coon dogs (Nyctereutes procyonoid.es) to humans (20). However, 
given that the human and civet SARS-CoV genomes were 96% 
identical and that the primary restriction to host range occurred at 
a small number of amino acid positions in the RBD, SARS-CoV 
may have been present in humans first and then been transmitted 
to civets and other animals in the markets (20). This is particularly 
intriguing, as the SARS-CoV takes a more generalist approach to 
host range, whereby it can infect cells expressing human, civet, 
and bat ACE2 (42, 43). In contrast, the civet SARS-CoV is re¬ 
stricted to infecting cells expressing civet ACE2 (42, 43). This 
study provides evidence that HCoV-NL63 is also a generalist, as it 
can infect bat cells as well as primate (48) and human (12) cells. 
Interestingly, HCoV-NL63 does not grow in mice. 

These are important observations because they indicate that 
some human and animal CoVs are generalists that can grow in a 
variety of different mammals, which is in contrast to many other 
mammalian CoVs that appear to be restricted to a single host. 
Recognizing that few mammalian coronaviruses have been tested 
in bat cell lines, these data suggest that humans and perhaps other 
select mammals may serve as the gateway species for the evolu¬ 
tionary expansion of mammalian CoVs. In addition, this observa¬ 
tion suggests that human and/or other mammalian CoVs could 
potentially circulate back and forth between bats and humans, 
establishing a zoonotic-reverse zoonotic cycle that may allow the 


virus to maintain viral populations in multiple hosts, exchange 
genetic information to alter pathogenesis or transmission charac¬ 
teristics, and potentially evolve variants that are capable of effi¬ 
ciently infecting humans or other mammals. Further work is nec¬ 
essary to determine the risks of such a cycle, evaluate the impact 
that this could have on public health, and determine whether or 
not this represents a general feature of the mammalian Corona- 
viridae, which are predicted to have originated from bats. 
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